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Abstract 

We calculate in closed form the complete 0(a%) color-singlet differential cross 
section for e + e~ — > 7* — > tpQ + X scattering. The cross section reduces at 
high energies to a heavy quark fragmentation form. We find that the energy 
scale at which the approximate fragmentation result becomes reliable exceeds 
the ipQ mass by more than an order of magnitude. We also discuss the color- 
singlet model's predictions for direct J/ip angular and energy distributions at 
CLEO. 
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I. INTRODUCTION 



During the past few years, there has been renewed interest in the study of heavy quarko- 
nium systems. Much of the recent work on this subject has been stimulated by large dis- 
crepancies between old predictions and new observations of Psi and Upsilon production at 
several experimental facilities. Orders of magnitude disagreements between theory and data 
have seriously undermined the conventional "color- singlet model" (CSM) picture of quarko- 
nia formation |T|-§[. In this model, charmonia and bottomonia mesons are presumed to 
exclusively originate from short distance processes that create heavy quark- ant iquark pairs 
in colorless configurations. The quantum numbers of pairs produced in high energy colli- 
sions on time scales short compared to Aqcd are required to precisely match those of the 
final state hadrons into which they nonperturbatively evolve. Although this CSM picture is 
simple, it does not explain several gross features of recent charmonia and bottomonia data 
collected at the Fermilab Tevatron |^§. It consequently must be abandoned as a complete 
theory. 

A new framework for treating quarkonia systems called Nonrelativistic Quantum Chro- 
modynamics (NRQCD) has been developed within the past few years 0. This effective 
field theory generalizes and improves upon the CSM in several regards. It allows for short 
distance processes to create heavy quark-antiquark pairs in color-octet configurations which 
can hadronize over much longer length scales into colorless final state quarkonia. Calcu- 
lations which include this color-octet mechanism appear to successfully describe Tevatron 
measurements |To|-[l3|. But in order to establish the validity of this new paradigm, it is 
necessary to consider quarkonia production in other experimental situations. 

Braaten and Chen have suggested that a clean signature of the color-octet mechanism 
may be observable in i/jq production at electron-positron colliders ||14j| . These authors have 
noted that the angular distribution of colored QQ pairs near the endpoint region may qual- 
itatively differ from that of their colorless counterparts. If this effect could be observed, 
it would support the color-octet production picture. It might also permit an independent 
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determination of the numerical values for certain NRQCD matrix elements. 

Before the search for color-octet quarkonia production in e + e~ annihilation can begin, 
one must first know the precise CSM prediction. Within the NRQCD framework, the color- 
singlet cross section is also expected to be quite accurate for all energies except near the 



endpoint region ]T^|. In this note, we therefore build upon previous studies reported in the 



literature |15|-|l8l and calculate the complete 0{a 2 s ) color-singlet cross section for e + e 



ipq + X scattering. We examine the contribution to ipQ production from the short distance 
modes e + e~ — > QQ[ 3 S[ lS) ]+g+g and e + e~ — > <5<3[ 3 ^i^]+Q + Q|] ari d we derive a closed form 
expression for the differential cross section. We then discuss the implications of the CSM 
result for direct J/ip observations at CLEO. Finally, we compare heavy quark fragmentation 
predictions with the color-singlet cross section and determine the energy scale at which 
fragmentation approximations become reliable. 

II. INCLUSIVE ANGULAR DISTRIBUTIONS IN ELECTRON-POSITRON 

COLLISIONS 

It is useful to note some general features of inclusive, unpolarized tpQ production in e + e~ 
annihilation. Unitarity, parity and angular momentum considerations restrict the form of 
the differential cross section expression 

AT f/ (e + (pi)e-( P2 ) - 7* - MPs) +X) = S(E 3 ) \l + a(E 3 ) cos 2 3 j . (2.1) 
dE 3 d cos 3 v / l j 

In particular, the allowed range for the angular coefficient function is constrained to lie 

within the interval —1 < a(E 3 ) < 1. We sketch a derivation of this result below. 

It is instructive to consider the subprocess J*(P) — > 4>q{pz) + X(P — p 3 ) where the 

intermediate photon is either longitudinally or transversely aligned. The squared amplitude 

for this decay 



1 We indicate the angular momentum and color-singlet quantum numbers of the QQ pair which 
hadronizes into the final state tpQ meson inside square brackets. 

3 



\A\ 2 = J2^(P^K(P;\TF^ (2.2) 

A 

involves a form factor F^ u which can be decomposed in terms of tensors that respect parity 
and gauge invariance: 

= -*i (<r - + J (rf - ^?^) (pa - ^p u ) ■ (2-3) 

Working in the 7* rest frame where the ipQ four-momentum looks like p 3 = (E 3 ,p 3 ) = 
(E 3 , \p 3 \ sin 9 cos 0, |p 3 | sin 9 sin 0, |p 3 | cos#), we find that the squared decay amplitude for a 
longitudinally polarized virtual photon reduces to 

\A\l = F 1 [l + a L ooa 2 0], (2.4a) 

with 

*L= lJ ^ry i - ( 2 - 4b ) 

For a transverse 7*, the squared amplitude takes the form 

\A\ 2 T = (2F 1 + ^F 2 ) [l+a T cos 2 #], (2.5a) 



where 



Since both |^4|| and |^4||. are nonnegative, eqs. (|2.4| ) and (|2.5|) imply > —1 and 
-1 < a T < 1. 

Helicity conservation requires the intermediate photon in e + e~ — > 7* — > ipQ + X to be 
transversely aligned relative to the beam axis in the m e = limit. The ipQ meson's angular 
distribution is therefore significantly restricted by simple symmetry considerations. In fact, 
the inclusive angular distribution of any unpolarized particle which is produced in electron- 
positron colliders operating well below the Z-pole goes as 1 + ct-r cos 2 9 with — 1 < a T < 1. 
So while observation of a pure sin 2 9 distribution for a lepton or hadron at a collider like 
CLEO is possible, a pure cos 2 9 distribution is not. As we shall see, all color-singlet ipQ 
predictions are consistent with these general considerations. 



III. COLOR-SINGLET i})q PRODUCTION 

The simplest parton level process which mediates color-singlet production of J PC = 1 
quarkonia is given by e + e _ — > QQ^S^} + g + g. Color, parity and charge conjugation 
conservation require two gluons to appear in the final state along with the colorless QQP'S'i^] 
pair. This channel consequently contributes to the ipQ cross section starting at 0{a 2 s ). Color 
singlet production also proceeds at the same order in perturbative QCD through the mode 
e + e~ — > QQ[ 3 S^] + Q + Q. These two distinct reactions have been considered separately in 
the literature |T5|-fr%|]. We will reexamine their joint impact upon Psi and Upsilon production 
and derive a closed form analytic expression for d 2 a/dE^ dcosO^. We can then compare 
the relative magnitudes of the gluon and quark processes as a function of center-of-mass 
energy 

The leading order diagrams which mediate e + (pi) e~{p2) — * QQ[ 3 S[^](p3) + g{pi) +g{p<o) 
and e + (pi)e~(p 2 ) — > QQ^S[ 1>> ](ps) + Q{p±) + Q{ph) scattering are illustrated in Figs. 1 
and 2. The hard collisions pictured in the figures form on short time scales a heavy quark 
and antiquark which fly out from the primary interaction point in nearly parallel directions 
and almost on-shell. The QQ pair then evolves over a much longer time interval into a 
physical i^q bound state. Working within the NRQCD framework and using computational 
methods discussed in Refs. P[ P^ , P^| , one can straightforwardly calculate the amplitudes for 



these processes. Their squares factorize into products of short distance coefficient functions 
and long distance NRQCD matrix elements. 

Integrating the squared amplitudes over the three particle phase space factor 

rf$ 3 = (2vr) + p 2 - P3 - p, - p 5 ) p (3- 1 ) 

is somewhat involved. As a simplifying measure, we rescale all dimensionful quantities rela- 
tive to the beam energy E and work with the dimensionless variables Z{ — Ei/E, = Pi/E, 
Xi = cos 9i and S = 2Mq/E. The phase space factor for the reaction with gluons in the final 
state can then be reduced to the form 
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QQ[ 3 s[ 1] ]( P3 ) 



e + {pi) 



e (Pi) 




g(P4) 



+ 



FIG. 1. Leading order Feynman graphs which mediate e + e — ► 7* 
scattering. 



(27r) 2 dzs dxs dz_ dw 
d<P 3 = — : — E 



where 



J{l-K*){\-xi)-w* 



permutations 



(3.2) 



Z- = 


Za — Z5 


<f-l = 


^4 - Az 3 + 6 2 + z 2 _ 


\Qs\ = 


v '4 - ^ 


K = 


Z-(2-z 3 ) 




w = 


X- + Kx 3 . 



(3.3a) 
(3.3b) 
(3.3c) 
(3.3d) 
(3.3e) 



The same result holds for the quark process with the simple alteration |<f_| = y 4 — 4^3 + z_. 

The available phase space volume clearly depends upon the masses of the final state 
bodies. For the e + e~ — > QQ^S^} + g + g channel, the limits of integration for the remaining 
energy and angular variables in eq. ( |3.2| ) are given by 

(3.4a) 
(3.4b) 
(3.4c) 

-,J(1 - K 2 ){1 -xl)<w< ^(1 - K 2 ){l-xl). (3.4d) 
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S< z 3 < 1 + - 
-1 < x 3 < 1 



-^4-5 2 <z_< Jzl-5 



\ 



"(pi) 
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+ 
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FIG. 2. Leading order Feynman graphs which mediate e + e — > 7* — 1 
— ► + X scattering. 

The corresponding limits for the e + e~ — > QQ[ 3 S[^] + Q + Q mode 

8 < z 3 < 1 
-1 < x 3 < 1 



QQ[ 3 S[ 1] ] + Q + Q 



l (4-4z 3 )(zl-Pj ^ < 



/ (4-4^ 3 2 -^)~ 
4 - 4^3 + <5 2 



1 - K 2 )(l - x\) < w < y/(l-K 2 )(l-x 2 3 ] 



(3.5a) 
(3.5b) 
(3.5c) 
(3.5d) 



are more tight due to the additional heavy quark and antiquark in the final state. 

After inserting the averaged squared amplitudes and reduced phase space factors into 
the formula 

1 



da = — -V \A\ 2 d<$> 



3 j 



(3.6) 



we can analytically integrate over w and Z- and obtain differential expressions of the form 
( |2.1|) . We display the resulting S(z 3 ) and a(z 3 ) functions for the e + e~ — > QQpSi 1 ^] +g + g 



ande+e" -> QQ[ 3 Sf 1} ] + Q + Q processes in the appendix. As a check, one can verify that 
|a g iuon| and |a quar k| do not exceed unity within their allowed z 3 ranges as required by the 
general constraints discussed in section II. The total 0(a 2 s ) angular coefficient function 

*Sgluon "gluon ~\~ Snu^-]^ Oquark / n \ 

"total = — — TT " (3-7) 

*J gluon \ >J quark 

also respects the bound — 1 < Qt to tai < 1- 

Another important check can be performed by considering the high energy behavior of 
the S(zs) and a(zs) functions. In the z% ^> 5 limit, the color-singlet cross section reduces to 

d 2 ° , + _ , , ^ 4tt (a s a EM Q Q ) 2 /„, ^ Q , 3l 



(eV^ Q U) = - ^ '7 (0|Or°( J gi)|0) l + cos^s (3.8) 
az-id cos 243 rrin^ x ' 1 J 



x 



z 3 (l - z 3 ) 2 (16 - 32z 3 + 72z 2 - 32zf + bz$) 



(2 - ^s) 6 

After integrating over cos 6 3 and recalling the relation (o|Cf Q ( 3 Si)|o) = 9|^(0)| 2 /2vr be- 
tween the color-singlet NRQCD matrix element and the i/jq wavefunction at the origin ||, 
we can write the ipq energy distribution as 

^(e+e- - ij Q + X) = 2a{e + e- -> QQ) x £> QW , Q (* 3 ) (3.9) 
az 3 



where Vq^ 1 p Q (z 3 ) denotes the heavy quark fragmentation function calculated in Ref. |H5 
The complete 0(a 2 ) color-singlet cross section thus correctly reproduces known fragmenta- 
tion results at high energies. 

IV. DIRECT J/ip PRODUCTION AT CLEO 



J/ip production is currently under study at CLEO pO| , pT| . Charmonia observed at this 
e + e _ facility mainly come from B meson decays. However, a clean sample of r/>'s originating 
from continuum production can be obtained by imposing a lower momentum cut on their 
dilepton decay products. Various characteristics of the resulting direct J/ip data sample 
can then be compared with predictions based upon color-singlet and color-octet production 



mechanisms. Such experimental investigations are underway |22 
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The angular distribution of direct J /if) mesons represents one observable which can be 
measured at CLEO. In Fig. 3, we plot the CSM prediction for the angular coefficient func- 
tion a. The results displayed in the figure are based upon the input parameter values 
E = 5.29 GeV, m c = 1.48 GeV, a s (2m c ) = 0.28, a EM (2m c ) = 0.0075, Q c = 2/3 and 
(0\O(^( 3 Si)|0^ = 1.2 GeV 3 . The dashed curve illustrates the function a g i uon associated 
with e + e~ — > cc[ 3 S , J 1 ' ) ] + g + g scattering. The shape of this curve agrees with numerical 

1 | I I I | I I I | I I I | I I I J I I n 



0.8 : 
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0.4 ; 

0.2 : 
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0.9 



1 



FIG. 3. Angular coefficient functions a g i UO n (dashed line), a quar k (dotted line) and a tota i (solid 
line) plotted against dimensionless energy variable 2:3. 



results of Driesen et al. reported in Ref . 0] . The dotted line in Fig. 3 depicts the function 
«quark originating from the e + e~ — > cc[ 3 S^] + c+c mode. The shape of a quar k is clearly quite 
different than that of a g i UO n- But since S^uaxk is substantially smaller than S'giuon at CLEO 
energies, it has only a small impact upon the total color-singlet function atotai which is rep- 
resented by the solid curve in Fig. 3. It is important to note that a to tai is predicted within 
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the CSM to be negative at the largest allowed values for z 3 . On the other hand, color-octet 
effects may render a total positive in the endpoint region |TJ|]. The angular distribution of 
the most energetic J/ip's at CLEO can therefore provide a valuable test of the color-octet 
mechanism. 

The energy distribution of direct J/^'s is another quantity which can be used to probe 
theories of quarkonia production. In Fig. 4, we display the separate contributions to da/dz^ 



from the e + e — > cc^S^} + g + g and e + e — > cc^Si 1 '] + c + c channels along with the 
total CSM prediction. The sensitivity of this energy observable to the charm mode is more 

2 




FIG. 4. Contributions to da/dz^ from the gluon (dashed line) and quark (dotted line) modes 
plotted against 03. The CSM prediction for the total direct J/tp energy distribution is represented 
by the solid curve. 



pronounced than that of the angular coefficient function. The areas underneath the dashed, 
dotted and solid curves respectively equal 0.74 pb, 0.07 pb and 0.81 pb. The quark process 
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thus contributes at the 10% level to direct J/ip production at CLEO. 

The e + e~ — > cc[ 3 Si ] + c + c mode is significantly phase space suppressed compared 
to e+e- -> cc[ 3 5i 1} ] + # + g at CLEO energies. As a result, its impact upon charmonia 
observables is minor. However, it is interesting to examine the relative importance of these 
two color-singlet channels as a function of center-of-mass energy. We plot in Fig. 5 the modes' 
separate contributions to the integrated J/ ip cross section along with their sum versus 
2E. We also display the integral of the charm quark fragmentation approximation (|3.8| ). At 
low energies, the charm quark mode is negligible compared to its gluon counterpart. At larger 
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FIG. 5. Integrated cross sections for the gluon (dashed line) and charm quark (dotted line) 
modes plotted as a function of y/S. The sum of the two is shown by the solid curve. The 
approximate charm quark fragmentation cross section is depicted by the dot-dashed curve. 

values of y/S, it becomes relatively more important. Finally, in the charm fragmentation 
limit y/S 3> m c , the quark mode dominates. 
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As can be seen in Fig. 5, the charm quark fragmentation curve rapidly asymptotes to 
the e + e~ — > cc[ 3 5'j 1 ' ) ] + c + c cross section. But it is important to note that the crossover 
point at which the rates for the charm and gluon modes become equal occurs around y/S ~ 
50 GeV. Consequently, the fragmentation approximation does not accurately reflect the 
total color-singlet cross section until V~S exceeds 2m c by more than an order of magnitude. 
This result for J/if> production at lepton colliders is quite different than that for hadron 
accelerators. Previous investigations have found that fragmentation approximations are 
reasonably trustworthy for production of ^'s at the Tevatron with p± > 10 GeV Hl2| , |T3|| . 
The moral we thus draw from this study of charmonia at CLEO is that the validity of 
fragmentation predictions must be carefully checked on a case-by-case basis. 



ACKNOWLEDGMENTS 



We thank David Politzer for helpful discussions. The work of P.C. was supported in 
part by a DuBridge Fellowship and by the Department of Energy under Grant No. DE- 
FG03-92-ER40701. The work of A.L. was supported in part by the DOE under Grant No. 
DE-FG03-92-ER40701. 



12 



APPENDIX: 



We list here the color-singlet functions S and a which enter into the differential cross 
section ( |2.1|) at leading order in both the perturbative QCD and NRQCD velocity expan- 
sions. The contributions from the e + e~ — > QQ[ 3 Sj 1 ' ) ] + g + g and e + e~ — > QQ^S^] + Q + Q 
channels are separately displayed. 



QQ[ 3 S{ 1) ]+ g + g mode: 



S, 



gluon 



7T (a s a EM Q Q f 



0|C»f Q ( 3 5 1 )|0 



216 



5 E 5 (z 3 - 2) 2 {2z 3 - 5 2 ) 3 {zj - 5 2 ) 
-<5 2 (4 + <5 2 )(48 + 485 2 + 135 4 ) + 32<5 2 (4 + <5 2 )(4 + 35 2 )z 3 



+ 8(32 - 565 2 - 24<5 4 + 5 e )zl - 16(32 + 4<5 2 + 3<5 4 ) 



+ 112(4 + 8 2 )zl - 128z§ (2z 3 -5 2 



6 2 



+ 



5 4 (4 - <5 2 )(48 + 965 2 + 135 4 ) - 325 4 (28 - 35 2 - 35 4 )z 3 
+ 85 2 (16 - 405 2 - 275 4 + 5 6 )z 2 3 + 165 2 (56 + 145 2 - 35 4 )z 3 3 
- 16(4 - 5 2 ){A + 55 2 )z 4 \ {4z 3 - 4 - 5 2 ) In 



2^ - 5 2 + 2 Jz 2 



S 2 



2z 3 - 5 2 - 2 x /z 2 



^gluon (^3) 



7T (a s a BM Qq) 5 



O|0f Q ( 3 5!)|O 



216 
x I 4 



5 E 5 (z 3 - 2) 2 {2z 3 - 5 2 ) 3 {zl - 5 2 ) 
5 2 (64 + 805 2 + 765 4 + 75 6 ) - 965 4 (4 + 5 2 )z 3 

- 8(32 - 405 2 - 445 4 - 5 & )z 2 3 - 165 2 (28 + 35 2 )zl 

+ 16(20 + 75 2 )zl - 128^] ( 2 ^3 - 5 2 )^3 

5 4 (4 - 5 2 )(4 + 5 2 )(4 + 75 2 ) - 325 4 (1 - 5 2 )(4 + 35 2 )^ 3 



- 85 2 (16 + 405" + 575 4 + 5*)z 2 + 16<5 2 (8 + 585 2 + 35 4 )z? 



+ 16(16 - 325 2 - 5<T)z 4 (Az 3 - 4 - 5 2 ) In 



2z 3 -5 2 + 2Jz 2 - 5 2 



(Ala) 



(Alb) 



x 



2z 3 -5 2 - 2^4 - 5 2 J S gluon (z 3 ) 



e+e- -> QQ[ 3 g| } ] + g + g mode: 
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Q , , vr (a s a EM Q Q ) 2 (o|^(^i)|o) 

x < ^ ^A + P-^ p- [~ 326 ^ + ^)( 48 + 2252 + 
+ 325 4 (768 + 400<5 2 + 665 4 + 35 6 )z 3 

- 165 2 (384 + 19205 2 + 556<5 4 + 295 6 - 25 8 )z 2 3 
+ 8<5 2 (1792 + 128<5 2 - 5685 4 - 805 6 - 5 8 )z 3 3 
+ 2(2048 - 110085 2 + 107525 4 + 31765 6 + 985 8 + 3<5 10 )4 

- 4(4096 - 78085 2 + 3424<5 4 + 6005 6 + 175 8 )z 5 3 
+ (38912 - 206085 2 + 4544<5 4 + 508<5 6 - 35 8 )z 6 3 

- 4(13312 - 8005 2 + 1205 4 - 35 6 )4 + 8(4512 - 205 2 - 155 4 )4 

- 32(336 - 5 2 )zl + 12804°] ( A2a ) 
8<5 4 (48 + 22<5 2 + 35 4 ) - 325 4 (24 + 55 2 )^ 3 

- 2<5 2 (448 + 165 2 + 85 4 - 35 6 )z 2 3 + 165 2 (56 - 105 2 - 55 4 )z| 
+ <5 2 (1152 + 2725 2 - 35 A )zl + 8(32 - 925 2 + 55 A )z 5 3 

- 56(16 + 5 2 )z % 3 +5124 



/ \ KP^S ®EM 

ttquarkl^3j — 



23 V4 + ^-4^3-2 v /(l-^)(4-^ 2 ) 
Qq)2 (01^(^)10) 



3888 <5 3 £ 5 z| (z 3 - 2) 6 (4 - 5 2 ) 

x < ^ ^A + P-^ P- [ 32 ^ 4(4 + 52)(16 + 2 ^ + 3 ^ 4) 

- 325 4 (256 + 485 2 + 225 4 + 35 e )z 3 

+ 165 2 (1152 + 10245 2 - 1405 4 - 535 6 - 25 8 )z 2 3 

- 85 2 (5376 + 1285 2 - 15765 4 - 2405 6 - 5 8 )z 3 3 
+ 2(2048 - 7685 2 - 199685 4 - 69685 6 - 3505 s - 35 10 )4 

- 4(4096 - 200965 2 - 111685 4 - 12085 6 - 43<5 8 )4 
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+ (38912 - 753925 2 - 169605 4 - 9965 6 - 3<5 8 )^ 

- 4(13312 - 63045 2 - 8725 4 - 35 6 )z 7 3 + 8(4512 - 5005 2 - 155 4 )^ 

- 32(336 - 5 2 )zl + 12804°] (A2b) 
+ [85 4 (16 + 25 2 + 35 4 ) - 325 4 (8 - 5 2 )z 3 

- 25 2 (320 - 2725 2 + 645 4 - 3<5 6 )z 2 + 165 2 (40 - 545 2 - 5<5 4 )zf 

- (1024 - 720<5 4 - 35 e )zl + 8(96 - 365 2 - 55 4 )z 5 3 
+ 8(80 + 75 2 )zl - 5124 

«*(4 - q'h * ^Z^l + JESSES ) x ' 



z 3 V4 + 5 2 - 4z 3 - 2^(1 -z 3 ) (4 -5 2 ) j 5 quark (z 3 ) 
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